A 3D model for carbon monoxide molecular line emission as a potential cosmic microwave background polarization contaminant

We present a model for simulating Carbon Monoxide (CO) rotational line emission in molecular clouds, taking account of their 3D spatial distribution in galaxies with different geometrical properties. The model implemented is based on recent results in the literature and has been designed for performing Monte-Carlo simulations of this emission. We compare the simulations produced with this model and calibrate them, both on the map level and on the power spectrum level, using the second release of data from the Planck satellite for the Galactic plane, where the signal-to-noise ratio is highest. We use the calibrated model to extrapolate the CO power spectrum at low Galactic latitudes where no high sensitivity observations are available yet. We then forecast the level of unresolved polarized emission from CO molecular clouds which could contaminate the power spectrum of Cosmic Microwave Background (CMB) polarization B-modes away from the Galactic plane. Assuming realistic levels of the polarization fraction, we show that the level of contamination is equivalent to a cosmological signal with r 720.02. The Monte-Carlo MOlecular Line Emission (MCMole3D) Python package, which implements this model, is being made publicly available

Characterization of foreground emission on degree angular scales for CMB B-mode observations: Thermal dust and synchrotron signal from Planck and WMAP data

We quantify the contamination from polarized diffuse Galactic synchrotron and thermal dust emissions to the B modes of the cosmic microwave background (CMB) anisotropies on the degree angular scale, using data from the Planck and Wilkinson Microwave Anisotropy Probe (WMAP) satellites. We compute power spectra of foreground polarized emissions in 352 circular sky patches located at Galactic latitude | b | > 20\ub0, each of which covers about 1.5% of the sky. We make use of the spectral properties derived from Planck and WMAP data to extrapolate, in frequency, the amplitude of synchrotron and thermal dust B-mode spectra in the multipole bin centered at \u2113 43 80. In this way we estimate the amplitude and frequency of the foreground minimum for each analyzed region. We detect both dust and synchrotron signal on degree angular scales and at a 3\u3c3 confidence level in 28 regions. Here the minimum of the foreground emission is found at frequencies between 60 and 100 GHz with an amplitude expressed in terms of the equivalent tensor-to-scalar ratio, rFG,min, between 3c0.06 and 3c1. Some of these regions are located at high Galactic latitudes in areas close to the ones that are being observed by suborbital experiments. In all the other sky patches where synchrotron or dust B modes are not detectable with the required confidence, we put upper limits on the minimum foreground contamination and find values of rFG,min between 3c0.05 and 3c1.5 in the frequency range 60-90 GHz. Our results indicate that, with the current sensitivity at low frequency, it is not possible to exclude the presence of synchrotron contamination to CMB cosmological B modes at the level requested to measure a gravitational waves signal with r 43 0.01 at frequency 72 100 GHz anywhere. Therefore, more accurate data are essential in order to better characterize the synchrotron polarized component and, eventually, to remove its contamination to CMB signal through foreground cleaning. \ua9 2016 ESO

Smoke in the Pipe Nebula: dust emission and grain growth in the starless core FeSt 1-457

J. Forbrich, et al., “Smoke in the Pipe Nebula: dust emission and grain growth in the starless core FeSt 1-457”, Astronomy & Astrophysics, Vol 580, August 2015. This version of record is available online at: https://doi.org/10.1051/0004-6361/201425375 Reproduced with Permission from Astronomy and Astrophysics, © ESO 2016.(abridged) Methods: We derive maps of submillimeter dust optical depth and effective dust temperature from Herschel data that were calibrated against Planck. After calibration, we then fit a modified blackbody to the long-wavelength Herschel data, using the Planck-derived dust opacity spectral index beta, derived on scales of 30' (or ~1 pc). We use this model to make predictions of the submillimeter flux density at 850 micron, and we compare these in turn with APEX-Laboca observations. Results: A comparison of the submillimeter dust optical depth and near-infrared extinction data reveals evidence for an increased submillimeter dust opacity at high column densities, interpreted as an indication of grain growth in the inner parts of the core. Additionally, a comparison of the Herschel dust model and the Laboca data reveals that the frequency dependence of the submillimeter opacity, described by the spectral index beta, does not change. A single beta that is only slightly different from the Planck-derived value is sufficient to describe the data, beta=1.53+/-0.07. We apply a similar analysis to Barnard 68, a core with significantly lower column densities than FeSt 1-457, and we do not find evidence for grain growth but also a single beta. Conclusions: While we find evidence for grain growth from the dust opacity in FeSt 1-457, we find no evidence for significant variations in the dust opacity spectral index beta on scales 0.02x36x30'). The correction to the Planck-derived dust beta that we find in both cases is on the order of the measurement error, not including any systematic errors, and it would thus be reasonable to directly apply the dust beta from the Planck all-sky dust model. As a corollary, reliable effective temperature maps can be derived which would be otherwise affected by beta variations.Peer reviewe

Planck 2015 results. V. LFI calibration

We present a description of the pipeline used to calibrate the Planck Low Frequency Instrument (LFI) timelines into thermodynamic temperatures for the Planck 2015 data release, covering four years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin-synchronous modulation of the cosmic microwave background dipole, but we now use the orbital component, rather than adopting the Wilkinson Microwave Anisotropy Probe (WMAP) solar dipole. This allows our 2015 LFI analysis to provide an independent Solar dipole estimate, which is in excellent agreement with that of HFI and within 1σ (0.3% in amplitude) of the WMAP value. This 0.3% shift in the peak-to-peak dipole temperature from WMAP and a general overhaul of the iterative calibration code increases the overall level of the LFI maps by 0.45% (30 GHz), 0.64% (44 GHz), and 0.82% (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power spectrum measured by WMAP. We estimate that the LFI calibration uncertainty is now at the level of 0.20% for the 70 GHz map, 0.26% for the 44 GHz map, and 0.35% for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration since the previous data release

Planck 2015 results. VI. LFI mapmaking

This paper describes the mapmaking procedure applied to Planck Low Frequency Instrument (LFI) data. The mapmaking step takes as input the calibrated timelines and pointing information. The main products are sky maps of I, Q, and U Stokes components. For the first time, we present polarization maps at LFI frequencies. The mapmaking algorithm is based on a destriping technique, which is enhanced with a noise prior. The Galactic region is masked to reduce errors arising from bandpass mismatch and high signal gradients. We apply horn-uniform radiometer weights to reduce the effects of beam-shape mismatch. The algorithm is the same as used for the 2013 release, apart from small changes in parameter settings. We validate the procedure through simulations. Special emphasis is put on the control of systematics, which is particularly important for accurate polarization analysis. We also produce low-resolution versions of the maps and corresponding noise covariance matrices. These serve as input in later analysis steps and parameter estimation. The noise covariance matrices are validated through noise Monte Carlo simulations. The residual noise in the map products is characterized through analysis of half-ring maps, noise covariance matrices, and simulations

Planck 2013 results. IX. HFI spectral response

The Planck High Frequency Instrument (HFI) spectral response was determined through a series of ground based tests conducted with the HFI focal plane in a cryogenic environment prior to launch. The main goal of the spectral transmission tests was to measure the relative spectral response (including out-of-band signal rejection) of all HFI detectors. This was determined by measuring the output of a continuously scanned Fourier transform spectrometer coupled with all HFI detectors. As there is no on-board spectrometer within HFI, the ground-based spectral response experiments provide the definitive data set for the relative spectral calibration of the HFI. The spectral response of the HFI is used in Planck data analysis and component separation, this includes extraction of CO emission observed within Planck bands, dust emission, Sunyaev-Zeldovich sources, and intensity to polarization leakage. The HFI spectral response data have also been used to provide unit conversion and colour correction analysis tools. Verifications of the HFI spectral response data are provided through comparisons with photometric HFI flight data. This validation includes use of HFI zodiacal emission observations to demonstrate out-of-band spectral signal rejection better than 10^8. The accuracy of the HFI relative spectral response data is verified through comparison with complementary flight-data based unit conversion coefficients and colour correction coefficients. These coefficients include those based upon HFI observations of CO, dust, and Sunyaev-Zeldovich emission. General agreement is observed between the ground-based spectral characterization of HFI and corresponding in-flight observations, within the quoted uncertainty of each; explanations are provided for any discrepancies.Comment: 27 pages, 28 figures, one of the papers associated with the 2013 Planck data releas

Planck 2015 results: XXV. Diffuse low-frequency Galactic foregrounds

We discuss the Galactic foreground emission between 20 and 100 GHz based on observations by Planck and WMAP. The total intensity in this part of the spectrum is dominated by free-free and spinning dust emission, whereas the polarized intensity is dominated by synchrotron emission. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with radio recombination line templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude H\u3b1 emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (\ue2\u2030 30%) of H\u3b1 having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak (in I\u3bd) ranging from below 20 GHz to more than 50 GHz. There is a strong tendency for the spinning dust component near many prominent H ii regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photo-dissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the Commander solution finds more anomalous microwave emission (AME) than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys, particularly at 5-20 GHz, will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck and WMAP data to make the highest signal-to-noise ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure. We argue that nearly all the emission at 40deg > l >-90deg is part of the Loop I structure, and show that the emission extends much further in to the southern Galactic hemisphere than previously recognised, giving Loop I an ovoid rather than circular outline. However, it does not continue as far as the "Fermi bubble/microwave haze", making it less probable that these are part of the same structure. We identify a number of new faint features in the polarized sky, including a dearth of polarized synchrotron emission directly correlated with a narrow, roughly 20deg long filament seen in H\u3b1 at high Galactic latitude. Finally, we look for evidence of polarized AME, however many AME regions are significantly contaminated by polarized synchrotron emission, and we find a 2\u3c3 upper limit of 1.6% in the Perseus region

Planck 2015 results. IV. Low Frequency Instrument beams and window functions

This paper presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to –25 dB at 30 and 44 GHz, and down to –30 dB at 70 GHz. It has been confirmed that the agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are 0.7% and 1% at 30 and 44 GHz, respectively (at ℓ ≈ 600); and 0.5% at 70 GHz (at ℓ ≈ 1000)

Planck 2013 results. XXXI. Consistency of the Planck data

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70 GHz (amplifier) and 100 GHz (bolometer) channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143 GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70, 100, and 143 GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8 to 15' are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement (measured by deviation of the ratio from unity) between 70 and 100 GHz power spectra averaged over 70 64 l 64390 at the 0.8% level, and agreement between 143 and 100 GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100, 143, and 217 GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematicerrors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described. Subsequent analysis of the beams quantifies the importance of asymmetry in the near sidelobes, which was not fully accounted for initially, affecting the 70/100 ratio. Correcting for this, the 70, 100, and 143 GHz power spectra agree to 0.4% over the first two acoustic peaks. The likelihood analysis that produced the 2013 cosmological parameters incorporated uncertainties larger than this. We show explicitly that correction of the missing near sidelobe power in the HFI channels would result in shifts in the posterior distributions of parameters of less than 0.3\u3c3 except for As, the amplitude of the primordial curvature perturbations at 0.05 Mpc-1, which changes by about 1\u3c3. We extend these comparisons to include the sky maps from the complete nine-year mission of the Wilkinson Microwave Anisotropy Probe (WMAP), and find a roughly 2% difference between the Planck and WMAP power spectra in the region of the first acoustic peak. \ua9 ESO, 2014